Quantum Spin Holography with Surface State Electrons 
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In a recent paper Moon and coworkers [C.R. Moon et a/., Nature Nanotechnology 4, 167 (2009)] have shown 
that the single-atom limit for information storage density can be overcome by using the coherence of electrons 
in a two-dimensional electron gas to produce quantum holograms comprised of individually manipulated 
molecules projecting an electronic pattern onto a portion of a surface. We propose to further extend the 
concept by introducing quantum spin holography - a version of quantum holographic encoding allowing to 
store the information in two spin channels independently. 
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In the era of ever increasing information storage den- 
sity the physical entity setting the limit to the progress 
has quickly been reduced to the size of a single atom. 
The concept of representing a single bit of information 
by a single atom or molecule has been conceived with 
the ability to manipulate individual atoms on a surface 
with the tip of a scanning tunneling microscope^ Two 
decades later it has also been shown that it is possible 
to store information in the state of a single free-standing 
atom trapped in an optical cavity 2 and that constructing 
truly atomic logic elements is not any more beyond the 
grasp of state-of-the-art experimental techniques 3 And 
despite the fact that achieving information storage den- 
sities close to the single-atom limit in real-world devices 
is still far from being feasible, the search for an informa- 
tion carrier which would break the limit set by the finite 
spacing between single atoms in solid state systems has 
not been abandoned. 

The main requirement that such a carrier has to ful- 
fill is that it's characteristic lateral dimensions should 
be smaller than those of a single atomic unit. In a re- 
cent paper Moon and coworkers have proposed that in- 
formation can be stored in a fermionic state of a coher- 
ent two-dimensional electron gas. 4 They have dubbed the 
proposed concept quantum holographic encoding, drawing 
parallels between the free space optical waves and elec- 
tron wavefuctions. They use atomic manipulation 1 to 
construct molecular holograms of carbon monoxide (CO) 
molecules on the (111) surface of a Cu crystal (Fig. [T|) 
hosting a two-dimensional quasi-free electron gas - a 
Shockley type surface state (SS)£ Electron waves scatter 
at local potentials introduced by the molecules and inter- 
fere (owing to inherent long range coherence 6 ) in a des- 
ignated area of the surface (white area in Fig. [1]) to form 
an electron density of states (LDOS) pattern represent- 
ing an information page. 4 This page can then be read out 
with a density sensitive scanning probe technique, such 
as the scanning tunneling microscopy (STM). There are 
several major advantages of such information encoding 



procedure: the information is projected onto a medium 
(i) free of lattice constraints and (ii) native to the system 
itself. This means that the size of a single bit is limited 
neither by the host surface lattice nor by the limitations 
of a readout technique (f.e. by lensing or collimation of 
an external readout bearn^). 

Furthermore, owing to the parabolic dispersion of the 
SS, the interference patters formed by electrons with dif- 
ferent energies are different. Using this property Moon 
et al. were able to project the hologram not only in two 
spatial degrees of freedom but also in the energy dimen- 
sion (Fig. [J), a concept very similar to optical volumetric 
holography* 4 ^ 

It is, however, conceivable, that in the same way as 
Moon and coworkers 4 have utilized the energy disper- 
sion, other intrinsic properties of electrons can act as a 
new dimension for information storage, pushing further 
back the information density limit. In the present let- 
ter we propose to use the spin of electrons as such ad- 
ditional dimension. We show that if the molecules or 
atoms used in the construction of a molecular hologram 
have pronounced magnetic properties or the surface state 
is inherently spin-polarized the scattering of surface state 
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FIG. 1. (color online) Principles of quantum holographic en- 
coding in a two-dimensional electron gas. Inspired by Ref. 0. 
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electrons becomes spin-dependent, allowing one to store 
different information pages in different spin channels of 
a collinear uniaxial system. As an example we discuss 
Cu(lll) surface decorated with bilayer cobalt (Co) is- 
lands capped with Cu thus obtaining a system with spin- 
polarized surface states. 8 With its help we demonstrate 
the possibility of simultaneously encoding different in- 
formation pages with electrons of the same energy but 
opposite spins as sketched in Fig. [TJ 

In order to have an effective means to treat the quan- 
tum holographic encoding theoretically we adopt the 
multiple scattering formalism, which has in the past 
proven to be a uniquely suitable tool for the descrip- 
tion of surface state scattering by impurities. 9 In terms 
of complex amplitudes the process of scattering at a sin- 
gle potential modifies the amplitude of electrons as 
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where ai n is the incoming amplitude and r] is the complex 
phase shift. The latter can be decomposed into real and 
imaginary parts such that e 2n? = aexp(2i(5), where the 
real part S describes the scattering phase and a describes 
the attenuation of the amplitude due to various processes 
(e.g. adsorption or, in particular case of an atom on a sur- 
face, scattering into bulk states). Parameters a and S are 
determined by the shape of the scattering potential and, 
in turn, uniquely define the process of scattering. The 
propagation of an electron across the surface can, to the 
first order, be considered as a free electron propagation 9 : 
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where k is the wave length and Ar is the propagation 
distance. The local density of states (LDOS) at a certain 
point r induced by a single atom located at is then de- 
termined by the interference of outgoing (in the spirit of 
the Huygens principle) and scattered amplitudes at that 
particular point. Considering the outgoing amplitude to 
be an order zero bessel function (being equal to 1 at the 
origin), the LDOS is then proportional to: 



LDOS(r,k) oc Re [a] = Re 
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The generalization to the multiple scattering case of n 
potentials is straightforward 9 -: 



a M s(r,k) = a Q [l - A] 1 a s 
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where a Q and a s are vectors of length n describing the 
propagation of the outgoing and scattered waves to and 
from scatterers i = l..n, and A is an n x n matrix describ- 
ing electron propagation between individual scatterers. 
The spatial LDOS distribution produced by an array of 
scatterers in a surface state is thus determined by 3 pa- 
rameters: the /c-vector of electrons (tied to their energy 
by the dispersion relation) and the parameters a and 6 



describing the scattering potentials. Moon et al have 
made use of the energy dispersion of the surface state 
(k = k(E)) to open the energy dimension for quantum 
holographic encoding. 4 We, however would like to uti- 
lize the fact, that some or all of the three above named 
parameters might be spin dependent. 

Two possibilities to achieve such spin-dependence 
spring to mind. One would be to use magnetic atoms or 
molecules for the construction of the hologram. The scat- 
tering potentials that they present for the surface state 
electrons shall then be different for electrons with differ- 
ently oriented spins. In a uniaxial system one could say 
that a = a(cr) and 5 = 5(a) become dependent on the 
spin a of the electron (either ^ or ^ in a uniaxial system), 
so that the same arrangement of scatterers would produce 
different LDOS patterns in different spin channels. In re- 
alization of such a scenario care should be taken to ensure 
the uniaxial character of the spins' alignment, their ther- 
mal stability and the absence of Kondo screening of the 
atomic or molecular moments. The first two and to some 
extent the third criterion could be satisfied, f.e., if one 
uses a substrate with substantial crystalline anisotropy 
and/or couples the spins of the adatoms or molecules 
constituting the hologram to an external magnetic field 10 
or an underlying subsurface magnetic layer iii2- To avoid 
the Kondo screening one could also concentrate on a sys- 
tem with a Kondo temperature below the experimental 
conditions. We, however, concentrate on another possi- 
bility of making electron scattering spin-dependent. In- 
stead of using magnetic statterers, we switch to a sys- 
tem with an inherently spin-polarized surface state - Co 
nanoislands or multilayers on Cu(lll). 5 This system is 
well studied both experimentally and theoretically^^ 
and provides us exactly with what we need - the spin- 
dependence of scattering parameters a and 5 as well as 
the spin-dependence of the electron wave vector k. How- 
ever, Co nanoislands do themselves scatter the surface 
state electrons^ and would thus slightly interfere with the 
encoding of the hologram. To avoid that, and at the same 
time protect the system from additional contamination 
and intermixing, the Cu surface with islands grown on 
top can be covered with a few additional layers of copper 
thus sealing the island within the surface and reducing 
their effect on the scattering of surface state electrons 
while still retaining the spin-polarized character of the 
surface state. Without limiting the generality we shall, 
in the following, concentrate on a particular example of 
such a system - bilayer cobalt nanoislands adsorbed on 
a Cu(lll) surface and capped with two additional lay- 
ers of copper [Cu/Co/Cu(lll)]. A sketch of the system is 
given in Fig.[2ja). We treat the system theoretically with 
a first-principles code based on Korringa-Kohn-Rostoker 
(KKR) Green's function method (in the framework of 
the density functional theory). This method is described 
in detail in numerous publications ^ 14 i 15 and has been 
shown to be well suited for describing the kind of sys- 
tems we are about to address £ 

The calculated spectral density maps of electrons some 
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FIG. 2. (color) (a) Sketch of the system: bilayer Co island on Cu(lll) covered with two monolayers of Cu. (b) Spectral density 
map of majority (left) and minority (right) electrons (along the K — T — K direction of the surface Brillouin zone) some 5 A 
above the system shown in (a), (c) Spin resolved Friedel oscillations at 0.6 eV around a singe Cu adatom on Cu/Co/Cu(lll). 
The data points represent the KKR calculation, the solid line is a fit of Eq. [3] to the data (modified for minority electrons to 
account for a faster decay). 



5 A above the surface of the system described above is 
presented in Fig. [2^b) for majority (left) and minority 
(right) electrons. It is apparent that the surface band 
structure of the system is spin-polarized. While majority 
electrons show a parabolic dispersion band (white dash- 
dotted line in the left panel) residing in the projected 
copper bulk band gap (black dashed lines), the minor- 
ity electrons (right panel) form a dispersive band below 
the Fermi energy with a relatively high effective mass 
and have a number of additional dispersive features in 
the energy ranges of 0.2 — 0.6 eV (arising from the hy- 
bridization of surface sp electrons with a d-level of Co) 
and above 0.7 eV (the analogue of the majority parabolic 
surface state band pushed beyond the Fermi level). It 
is worth noting, that those bands cross the boundary 
of the projected bulk band gap, which in normal cir- 
cumstances would lead to a formation of a fairly broad 
resonance. Here, however, the bands remain largely un- 
broadened, which means that the layered structure topo- 
logically protects the surface state from scattering into 
the bulk. This protection is, however, not complete as 
shall be seen later. 

Having at our disposal a spin-polarized surface state we 
can now attempt to encode a hologram using electrons 
of a certain energy but with different spin orientations to 
store different information pages. To do so, we need to 
extract the scattering parameters (a, S and k). In the 
selection of the electron energy for holographic encoding 
we rely on two factors: high k- vector to reduce the char- 
acteristic electron wavelength and thus maximize the po- 
tential information density and possibly small deviation 
from the Fermi energy in order to facilitate experimental 
realization. Based on that we select an energy of 0.6 eV 
above the Fermi energy, which corresponds to the upper 
end of the dispersive minority band in Fig.[2fb). Param- 
eters a, S and k are then obtained for the given energy by 
fitting Eq.[3]to a radial LDOS distribution (Friedel oscil- 



lations) produced by a single scatterer on a surface. Such 
a distribution for a single Cu adatom on Cu/Co/Cu(lll) 
is given in Fig. [2fc) for majority (red circles) and minor- 
ity (blue rectangles) electrons. It is apparent that the 
Friedel oscillations of majority and minority electrons are 
substantially different. Not only are they denned by dif- 
ferent scattering parameters, but they also seem to have 
different decay ratios. While the majority curve can be 
well fitted with Eq. yielding the red solid line, the mi- 
nority curve seems to decay exponentially rather than 
follow the 1/r rule 16 . This is the consequence of the up- 
per part of the minority band (at 0.6 eV), with which the 
scattered electrons can be associated, lying outside the 
projected copper bulk band gap. While the Co bilayer 
island and the capping Cu bilayer seem to largely pre- 
vent the surface state electrons from being scattered into 
the bulk, the probability of such a scattering is non-zero 
and thus the coherence length of surface state electrons is 
reduced causing the exponential, rather than 1/r, decay 
of Friedel oscillations. The blue solid curve represents 
a fit of the data to Eq. [3] with an additional exponen- 
tial factor added. Nevertheless ,the conditions for the 
spin-selective encoding are satisfied and we can assume 
that encoding different information pages with electrons 
of different spin character is possible. Following the no- 
tation of Moon and coworkers such encoding could be 
called the "quantum spin holography" (QSH). 

To prove the possibility of QSH we will try to encode 
the patterns used by Moon et al. in the "volumetric 
quantum holography" section of their paper 4 - namely 
5 x 7 bit pages in likeness of letters "S" and "U". The 
procedure of encoding the quantum spin hologram is vir- 
tually identical to the one described by Moon et al for 
the volumetric quantum holography: one has to find a 
scatterer distribution producing the desired pattern in 
the target area of the surface (Fig. [1]). The pattern pro- 
duced by a given atom distribution (atomic hologram) 
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FIG. 3. (color) Quantum spin hologram, encoding two information pages in different spin channels. Bit maps 5x7 pixels in 
likeness of letters "S" (a) and "U" (b) encoded with majority and minority electrons, respectively. The insets show target (top) 
and resulting (bottom) information pages. On the right (c) the total LDOS distribution p(tot) — p(EfA) + p(Ef,-1) is shown. 



was calculated using the multiple scattering formalism 
as formulated in Eq. [4] The corresponding parameters 
a, S and k were obtained by fitting as described above. 
The quality of the pattern created by a particular distri- 
bution was assessed by calculating the standard Pearson 
correlation^ between the calculated LDOS pattern [dig- 
itized on a dense real space mesh and averaged locally to 
obtain the 5 x 7 bit pattern, see lower insets in Figs. [3] 
panels (a) and (b)] and the target bit map [upper insets 
in Fig. [3ja) and (b)]. The correlation was then maxi- 
mized over the whole set of possible atomic/molecular 
configurations to optimize the hologram. While theo- 
retical annealing was effectively utilized by Moon and 
coworkers^ to predict optimal molecular positions for vol- 
umetric quantum holography, we chose make use of a ge- 
netic algorithm technique - a search heuristic that mim- 
ics the process of natural evolutio n 18 ! 19 to enhance the 
theoretical annealing method. Genetic algorithms allow 
faster convergence and increase the chance of avoiding 
local minima in the optimization of the pattern-target 
correlation. In optimizing the atomic configuration of 
the molecular hologram, the following experimental ne- 
cessity was taken into account: the atoms were prohib- 
ited to form dimers as those are virtually impossible to 
control with an STM tip. 

The procedure described above applied to an area of 
Cu/Co/Cu(lll) surface approximately 10 x 10 nm with 
a target area of 2 x 3 nm yielded the holographic pat- 
terns shown in Fig. [3](a) for majority and in Fig. [3fb) 
for minority electrons at the Fermi level. The optimized 
positions of atoms forming the hologram can be deduced 
from the map of the total density of states in Fig. [3jc). 
It is apparent that while the total LDOS map [Fig. [3fc)] 
contains a seemingly random pattern of spots the ma- 
jority and minority [Figs. [3fa) and[3fb)] LDOS patterns 
closely resemble (within the projection area in the cen- 
ter) the targeted information pages (upper insets). To 



quantify the information contained in the hologram we 
apply the procedure used by Moon et al. and dividing 
the projection area into a grid of 5 x 7 squares average the 
LDOS within each square. The obtained values can then 
be compared to a predefined threshold resulting in 35 
bits of encoded information). The result of averaging of 
the holograms presented in panels (a) and (b) of Fig. [3] is 
presented in the lower insets in the corresponding panels. 

In summary, we have demonstrated the feasibility of 
quantum spin holography, i.e. information can be stored 
in the local density of states patterns created indepen- 
dently by majority and minority electrons scattered at 
an ensemble of magnetic atoms. Introducing magnetism 
into the system thus theoretically allows one to double 
the information density as compared to the volumetric 
quantum holographic encoding. 4 We are also rather op- 
timistic about the possibility of the experimental realiza- 
tion of the proposed concept using a spin-polarized scan- 
ning tunneling microscope for construction and readout 
of quantum spin holograms. 



1 B. Eigler and E. Schweizer, Nature 344, 524 (1990). 

2 H. P. Specht, C. Nolleke, A. Reiserer, M. Uphoff, E. Figueroa, 
S. Ritter, and G. Rempe, Nature 473, 190 (2011)] 

3 A. J Heinrich, C. P. Lutz, J. A. Gupta, and D. M. Ei- 
gler, Science New Series, 298, pp. 1381 (2002) ; A. A. Kha- 
jetoorians, J. Wiebe, B. Chilian, and R. Wiesendan- 
ger, |Science 332, 1062 (2011)| A. Heinrich and S. Loth, 
Science 33 2, 1039 (2011)1 

4 C. R. Moon, L. S. "Mattos, B. K. Fost er, G. Zeltzer, and 
H. C. Manoharan, |Nat Nano 4, 167 (2009)] E. J. Heller, 

JNat Nano 4, 141 (2009j] ~ 

°H. Oka, P. A. Ignatiev, S. Wedekind, G. Rodary, L. Nieber- 
gall, V. S. Stepanyuk, D. Sander, and J. Kirschner, 
Science 327, 843 (2010) ; N. N. Negulyaev, V. S. Stepanyuk, 
L. Niebergall, P. Bruno, M. Pivetta, M. Ternes, F. Patthey, 
and W.-D. Schneider, |Phys. Rev. Lett. 102, "24610 2 (2009) ; 
W. Shockley, |Phys. Rev. 56, 317 (1939j] ~ 

6 F. Silly, M Pivetta, M. Ternes, F. Patthey, J. P . Pelz, and W.-D. 
Schneider, |Phys. Rev. Lett. 92, 016101 (20047] 



5 



7 J. F. Heanue, M. C. Bashaw, and L. Hesselink, 
JScience 265, 749 (1994)] 

8 L. Diekhoner, M. A. Schneider, A. N. Baranov, V. S. Stepanyuk, 
P. Bruno, and K. Kern, |Phys. Rev. Lett. 90, 2368 01 (2003) 

9 E. J. Heller, M. F. Crommie, C. P. Lutz, and D. M. Eigler, Nature 
369, 464 (1994). 

1 F. Meier, L. Zhou, J. Wiebe, and R. Wiesendanger, 

|Science 320, 82 (2008)1 
11 0. O. Brovko, P. A. Ignatiev, V. S. Stepanyuk, and P. Bruno, 
JPhys. Rev. Lett. 1 01, 03680 972008)] 

V ^T. Uchihashi, J. Zhang, J . Kroger, and R. Berndt, 

|Phys. Rev. B 78, 033402 (2008) 
13 TT Okl K. T^c~ S7 



V. 



Stepanyuk, D. 



Wedekind, 
Sander, and 



G. Rodary, 
J. Kirschner, 



|Phys. Rev. Lett. 107, 187201 (2011)] 



14 K. Wildbe rger, V. S. St epanyuk, P. Lang , R. Zeller, and P. H. 

Dederichs, |Phys. Rev. Lett. 75, 509 (1995) | 
15 R. Zeller, P. H. D ederichs, B. U jfalussy, L. Szunyogh, 

and P. Weinberger, Phys. Rev. B 52, 8 807 (1995); N. N. 

Negulyaev, V. S. Stepanyuk, L. Niebergall, P. Bruno, 

W. Hergert, J. Rep p, K.-H. R ieder, and G. Meyer, 
JPhys. Rev. Lett. 1 01, 22660 1 (2008) | 

16 In fact, since the decay path cannot be unambiguously identified 
from our calculations, the correct analytic decay rate cannot be 
guaranteed. The exponential factor, however, yields much more 
accurate fitting results, than l/r n po wer law. 

17 J. L. Rodgers and W. A. Nic ewander, |Am. Stat. 42, 59 (1988)] 

18 S. V. Dudi y and A. Zunger,|Phys. Rev. Lett. 97, 046 401 (2006) 

19 S. Forrest, [Science 261, 872 (1993)| 



